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Abstract. The first purpose of this paper is to shed some new hght on the old 
question of selecting the number of beams in intensity-modulated radiation therapy 
(IMRT). The second purpose is to ilhiminatc the related issue of discrete static beam 
angles vs. rotational techniques, which has recently re-surfaced due to the advancement 
of volumetric arc therapy (VMAT). A specific objective is to find analytical expressions 
that allow one to address the points raised above. To make the problem mathematically 
tractable, it is assumed that the depth dose is flat and that the lateral dose profllc can 
be approximated by polynomials, specifically Chebyshev polynomials of the first kind, 
of finite degree. The application of methods known from image reconstruction then 
allows one to answer the first question above as follows: The required number of beams 
is determined by the maximum degree of the polynomials used in the approximation 
of the beam profiles, which is a measure of the dose variability. There is nothing to 
be gained by using more beams. In realistic cases, in which the variability of the 
lateral dose profile is restricted in several ways, the required number of beams is of 
the order of 10 to 20. The consequence of delivering the beams with a 'leaf sweep' 
technique during continuous rotation of the gantry, as in VMAT, is also derived in 
analytical form. The main effect is that the beams fan out, but the effect near the 
axis of rotation is small. This result can serve as a theoretical justification of VMAT. 
Overall the analytical derivations in this paper, albeit based on strong simplifications, 
provide new insights into, and a deeper understanding of, the beam angle problem in 
IMRT. The decomposition of the beam profiles into well-behaved and easily deliverable 
smooth functions, such as Chebyshev polynomials, could be of general interest in IMRT 
treatment planning. 
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1. Introduction 

The problem of selecting the number and direction of beams in intensity-modulated 
radiation therapy (IMRT) has been investigated by many and has been the subject of 
debate as early as in 1995 (Mohan & Ling 1995). One may argue that even today, 
after many more papers have been published on this topic, the results are not entirely 
convincing, and not practically useful overall. In fact, in today's clinical practice, the 
number and direction of beams often has to be found by trial and error. Two 'schools' of 
IMRT delivery co-cxist today: in the more common approach, a relatively small number 
(typically of the order of 10 or less) of fixed intensity-modulated beams is delivered with 
a multileaf collimator (MLC). In another approach, a very large number (> 50) of 
beams is used in the rotational tomotherapy approach (Mackie et al. 1993, Carol 1994). 
More recently, another rotational IMRT approach, 'Single-Arc', has been developed 
(Cameron 2005, Ulrich et al. 2007, Otto 2008, Wang et al. 2008) and has found a lot 
of interest as a commercial product (RapidArc, VMAT). In the Single- Arc approach 
there is no intensity modulation for any one beam angle, but the radiation field shape is 
varied dynamically and rapidly with an MLC as the gantry rotates around the patient. 
In the rotational approaches, the selection of beam angles is not a problem, but the 
distribution of dose over large vohimes of healthy tissues has been a concern. 

Intuitively one would expect that using more beams will always help to shape 
radiation dose distributions to match the tumor target volume, even though it is clear 
that there is a point of diminishing return. In this paper I will address the 'how many 
beams' problem from a theoretical viewpoint. I will introduce strong simplifications 
to make the problem mathematically tractable. The mathematical approach used 
here is similar to the one from an earlier paper on image reconstruction (Bortfeld 
& Oelfke 1999). One of the results is that, contrary to the intuition, there is no 
benefit whatsoever in increasing the number of beams beyond a certain threshold, which 
depends on the achievable amount of intensity modulation per beam. I will present some 
illustrative examples. 

In the last methodological chapter I will take the step from standard IMRT with 
a small number {K) of fixed beam angles towards rotational therapy. Prom what I 
said above it is clear that dynamically rotating the beam cannot yield better dose 
distributions, but it can potentially reduce the treatment time. First I will derive 
a simple equation to describe the effect of rotating the gantry while IMRT is being 
delivered with a sweep technique (Convery & Rosenbloom 1992). I will apply this 
equation to the delivery of polynomial beam profiles over short arcs of 180°/fC. As a 
result, the impact of this small angle rotation is small, thus serving as a theoretical 
justification of the Single- Arc approach. 
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2. Polynomial approximations 

Let us adopt the coordinate system from an earlier paper (Bortfeld & Oelfke 1999). 
Here the x-axis points from left to right, and the y axis points to the radiation source 
at the = 0° gantry position. In other words, at = 0° gantry angle, the beam goes in 
the —y direction, i.e., downwards. The gantry angle is counted positive in the counter- 
clockwise direction; note, however, that the gantry angle was defined in reversed (i.e., 
in clockwise) orientation in a more recent paper (Bortfeld & Webb 2009). 

I will use a zero-order approximation of the depth dose distribution, i.e., I assume 
that the depth-dose is fiat. This means in particular that I ignore beam divergence 
and the variability of scatter as a function of depth. Incidentally, one can consider 
this approximation as a first order approximation (linear depth-dose profile) if one uses 
parallel-opposed beams, but parallel-opposed beams will not be considered further in 
the rest of this paper. Within this approximation one can now write the two-dimensional 
dose resulting from a beam impinging under an angle with fiuence (intensity) profile 
Uip) as 

D4r) = Uif.n^). (1) 

Here r = (x, y) and n<^ = (cos 0, sin 0) is the unit vector perpendicular to the beam 
direction. In image reconstruction terminology, equation ([T]) is a 'backprojection', i.e., 
a 'smearing out' of fiuence values along the rays of the beam (Bortfeld & Oelfke 1999). 
In section [7] I will discuss why I believe that, in spite of this strong simplification, the 
results are not only of theoretical but also of some practical relevance. 

Let us now make one additional assumption, namely that one can write the fiuence 
(intensity) profiles as a set of polynomials of degree m. This should always be possible 
because in practice the fiuence profile is well behaved thanks to the finite source size, 
scatter, partial transmission through the leaf ends of the MLC, among other reasons. 

Under these assumptions, it follows from an earlier paper (Bortfeld & Oelfke 1999) 
that one can produce any dose distribution of the form 

M~l ^ 

D{r} =J2 ^M) (r- n^rd<P, (2) 

m=0 Q 

within the unit circle |r| < 1 using K > M evenly spaced beams. Here Wm{4>) 
are arbitrary weight factors. Equation ^ can be considered as a special form of a 
polynomial approximation of an arbitrary two-dimensional dose distribution D{f). A 
perhaps more important result of that paper is that the required number of beams equals 
M, the maximum degree polynomial plus 1. In particular, there is no benefit at all in 
using more than M beams. For a better understanding of these finding I will consider a 
specific set of polynomials, Chebyshev polynomials. 
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2.1. Chebyshev polynomials of the first kind 

Chebyshev polynomials of the first kind, T^ip) are a special set of (orthogonal) 
polynomials that can also be defined using cosine functions, namely: 

^m(p) = cos(m arccosp) (3) 

They are defined on the interval p G [—1,1] and can be considered as cosine or sine 
functions for even and odd m, respectively, with a frequency that increases away from the 
origin. On the other hand, they are also polynomials, for example T4(p) = 8p^ — + 1, 
Ts (p) = 128/ - 256/ + 160p^ - 32/ + 1 and Tg (p) = 256/ - 576p^ + 432^^ - 120^^ + Qp. 
Figure [l] shows those Chebyshev polynomials T4(/, Tg{p), and Tq{p). 




Figure 1: Examples of Chebyshev polynomials of the first kind, T^ip), for m = 4 (solid 
line), m = 8 (dotted line), and m = 9 (dashed line). 

In an earlier paper (Bortfeld & Oelfke 1999) some properties of Chebyshev 
polynomials of the second kind, Um{p), were utilized. Those Chebyshev polynomials of 
the second kind are Eigenf unctions of the operation consisting of a backprojection and 
subsequent projection under a different angle. The two kinds of Chebyshev polynomials 
are related through the following equation: 

TM = "-W-^-^W . (4) 

Here I will focus on Chebyshev polynomials of the first kind because they oscillate 
between —1 and +1, and issues with negativity (Cormack & Quinto 1989) can therefore 
be avoided simply by a zero-offset of +1. Chebyshev polynomials have been used in the 
theory of IMRT by Cormack et al. (Cormack 1987, Cormack & Cormack 1987, Cormack 
& Quinto 1989) but, to the best of my knowledge, they have not been used to get a 
handle on the beam angle problem. 
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3. 'Making' rotating beams from fixed beam angles 

The fact that one does not need more than M beams if each of the beam fluence profiles 
can be described by a polynomial of degree up to m = M — 1, will now be illustrated 
using the polynomial T/^{p) (the solid line in figure [T]) as an example. If one uses K = b 
fixed beams evenly spaced at angles (pk = 180° k/K, where k = 0, . . . , K — 1, and each 
of those 5 beams delivers a T4 fiuence profile of variable weight, then one can 'make' 
a beam with a T4 fiuence profile from an arbitrary angle (p. The weights of the fixed 
beams to be used are given by 

/ , N 1 sin((m + l)((/)fc-(/>)) 
K sm{(j)k - (j)) 

with m = 4 and K = 5 in this case. This result follows from appendix A in (Bortfeld & 
Oelfke 1999), in which Um has to be replaced with Um-2, and then from equation Q. 
If sin(0fc — 0) = 0, then wi^cpk) = ±(m + l)/K, and the sign can be determined with 
rHopital's rule. 

Figure |2] illustrates these findings graphically for = 0°, = 15°, and = 30°. The 
weights calculated with equation ^ are attached to the 5 beams in the plots. In figure 
[2^ with the beam in vertical direction (0 = 0°), only the vertical fixed beam (0o = 0°) 
has a weight of 1, and the four others a weight of 0. In figures [2)3 and[2]3 all five beams 
contribute with a finite weight to 'make' those oblique beam directions. 

For each of the beams, a zero-offset (i.e., a uniform fluence proflle) has to be added 
to make the fluence non-negative. The total zero-offset over all beams depends on the 
weights w{4>k), as calculated using equation ([s]), and therefore on 0. The biggest total 
zero-offset is required for = 180°/2i^' = 18°, and is close to 2. I therefore use a total 
zero-offset of 2 in all cases shown in flgure [2] I should mention that a T4 dose proflle with 
arbitrary orientation can also be generated from 5 flxed beams with intensity proflles 
that are polynomials of degree 4 but not Chebyshev T4 proflles. In other words, there 
is a 'nuUspace' that leads to non- uniqueness (Cormack 1987). However, one can show 
numerically that the use of 5 T4 proflles (evenly spaced) requires the smallest zero-offset. 

Moreover, it is interesting to note that beams with arbitrary orientation can also 
be generated from a number of flxed beams that are not evenly spaced but arbitrarily 
located. This flnding may serve as an explanation that the choice of beam angles 
is generally less critical in IMRT than in standard 3D conformal radiation therapy 
without intensity modulation. However, with non-equally spaced beams, the solutions 
(the equivalent of equation ([s])) may become highly oscillatory (i.e., one beam may 
require a large positive weight and the next a large negative weight), and very large 
zero-offsets may be required. 



4. A simple two-dimensional example 



For a further illustration of the flndings from above, let us now look at a simple 2-D 
IMRT example case: a circular target volume that 'wraps around' a circular critical 
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Figure 2: Making rotating beams from fixed beam angles. By varying tlie weiglits 
of intensity-modulated beams from fixed beam directions, one can create beams from 
arbitrary directions, (a): 0°, (b): 15°, (c): 30°. In every case, 5 beams from fixed 
evenly spaced directions (every 36°) are used, as shown by the white lines. Every beam 
produces a Chebyshev T4(p) fluence profile but with different weights that are calculated 
with equation ([s]) and shown as white numbers next to the beams. Negative values of 
the fluence profiles are avoided by a zero-offset, i.e. by adding constant fluence. See the 
text for further information. 
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structure. To match the region of support of the Chebyshev polynomials, the radius 
of the target is set to 1 unit, with its center positioned at {xo,yo) = (0,0). The 
critical structure has a radius of 1/3 units and its center is at (0, —2/3). This geometry 
corresponds with the geometry used in a recent paper (Bortfeld & Webb 2009). The 
dose prescription has a dose level of 1 unit to the target volume, and a reduction to 0.5 
units in the critical structure. The 50% dose reduction in the critical structure is chosen 
as a challenging but realistic value, and to avoid any issues with negativity. 

The calculation of the IMRT fiuence profiles is done with the recipe from an earlier 
paper, see page 1109 in (Bortfeld & Oelfke 1999). Here I approximate the fiuence profiles 
with M = 10 Chebyshev polynomials with degrees between m = and m = 9, see figure 
[TJ The specific steps of the calculation are: 

(i) Projection of the prescribed dose along K evenly spaced beam directions with angles 
0fc at lateral positions pj = cos(;^^|^) , j = 1, . . . , M. 

(ii) Discrete sine transform of each of those K projections. 

(iii) Multiplication of the coefficients with a 'ramp filter'. 

(iv) Composition of 'filtered projections' to obtain the K fiuence profiles. 

The last step above involves the use of Chebyshev polynomials of the second kind. For 
consistency with the rest of this paper, the algorithm can be re-written for the use of 
Chebyshev polynomials of the first kind through equation (|4]). 

After the calculation of the K fiuence profiles, those profiles are then backprojected, 
which simulates an irradiation with K beams. Figure [3]^a) shows the resulting dose 
distribution for K = 100 beams that are evenly spaced at every 1.8°. Now, as I said 
above, for polynomial fiuence profiles of degree up to m = 9, i^' = 10 beam angles should 
actually suffice. To illustrate this point, I use i^' = 10 beams at every 18° in figure [3]d. 
The 10 fiuence profiles are calculated from the case with 100 beams by use of equation 
(|5|. I should note that after the re-distribution into the 10 fiuence profiles, there is 
no guarantee that all of them are non-negative, even if the original 100 profiles are all 
non-negative. However, in this specific case negativity is not an issue. By comparing 
figures [3]^a) and|3](b) one can now see that the two dose distribution are indeed identical 
within the unit circle, which encompasses the target volume and the critical structure. 
Only outside of the unit circle the two distributions begin to differ more and more: with 
K = 100 beams the distribution is more uniform, whereas the case with i^' = 10 beams 
exhibits the typical ragged distribution. With fewer than M = 10 beams, there would 
be differences also within the unit circle. 

5. How many beams are needed? 

Can one now find a general answer to this question? As 1 said above, the variability of the 
fiuence profiles, measured by the degree of the Chebyshev polynomials approximating 
those fiuence profiles, determines how many beam are needed. The fiuence variability 
depends, in turn, on the complexity of the dose prescription, and is therefore case 
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10 beams 
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Figure 3: This figure shows dose distributions for the treatment of a circular target 
volume (white solid line) that 'wraps around' a small circular critical structure (white 
dashed line) with K = 100 beams (a) and K = 10 beams (b). The dose prescription 
is one dose unit to the target volume, and a dose reduction to 0.5 units is requested 
in the critical structure. All beams are composed of M = 10 Chebyshev polynomials 
with degrees between m = and m = 9. Within the target circle, which encompasses 
the critical structure, the two distributions in (a) and (b) are identical, i.e., there is no 
benefit in using more than 10 beams here. 
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dependent. However, even the most complex cases do not benefit from an arbitrarily 
high number of beams. This is so because there is also a physical limit on the achievable 
amount of fluence (or rather dose) modulation, which is due to the reasons mentioned 
above (scatter, transmission, finite source size). This physical limit reveals itself through 
the fact that the dose fall-off at a field edge is never infinitely steep, but exhibits a dose 
'penumbra'. Specifically, the width of the penumbra (80% to 20%) near the leaf-end of 
an MLC at a typical treatment depth is of the order of 7 mm (Boyer et al. 1992). As a 
consequence, the minimum width (full width at half maximum, FWHM) of a deliverable 
dose peak, Ap, is of the order of 1.4 x 7 mm = 1 cm (Bortfeld et al. 2000). Now, by 
interpreting Chebyshev fluence profiles as a series of dose peaks (see figure [T]), one can 
correlate this physical limit with the maximum degree polynomial that can be delivered, 
and thus with the required number of beams. 

Chebyshev polynomials Tmip) of the first kind have m roots. The roots cut the 
[—1,1] interval into m + 1 sub-intervals, which represent the peaks and valleys of Tm{p)- 
The average width of the sub-intervals is 2/(m + 1), hence the average peak width is 
also 2/(m + 1) (FWHM). To give this expression a quantitative meaning, one has to 
identify the 2 units in the numerator with 2R, where R is the radius of the the circular 
region of interest, which was set to 1 in the abstract examples above. The region of 
interest contains the target volume and the most important nearby critical structures, 
which may extend into the target volume. Thus the average peak width is 2R/ (m + 1). 
If one equates this with the physical limit, Ap above, and utilizes the fact that m + 1 
governs the required number of beams, K, one obtains: 

(6) 

Lip 

To give a couple of examples, for relatively small cases where the radius of the region 
of interest is i? = 5 cm, = 10 beams are required, whereas for larger cases with 
R = 10 cm, one needs K = 20 beams. For 'sharper' beams with smaller penumbra, i.e., 
with a smaller Ap, the number is correspondingly higher. 

I should emphasize that equation ^ can only serve as a rough estimate of the 
required number of beams, as will be discussed in section [7| 

6. Towards a theoretical justification of Single- Arc 

Now that it has been shown that a finite and moderate number, K, of fixed IMRT 
beams suffices, how can one justify Single- Arc IMRT (RapidArc, VMAT), in which the 
number of beams is effectively infinite, but without intensity modulation from any one 
beam direction? One way to provide a theoretical foundation of Single- Arc is as follows: 
Let us start with K beams and assume that the fiuence profile for each beam is delivered 



with a 'sweep' technique (Convery & Rosenbloom 1992), see Appendix A, in which the 
leaves move uni-directionally, say from left to right. What is the impact on the fiuence 
and dose distribution, if one delivers the dose from each beam not while the gantry is 
static at a fixed beam angle, but rotating over the angle 180° / K7 Some investigations to 
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answer this question were presented in (Webb & McQuaid 2009). A closed form answer 
is derived in Appendix A Basically the impact of gantry motion during leaf sweep IMRT 
delivery is a stretching or compression of the fiuence profil^ depending on whether the 
gantry motion increases (prograde motion) or decreases (retrograde motion) the motion 
of the field edges (the projected leaf ends). 

For leaf sweep IMRT, the fiuence profiles f{p) are typically decomposed into two 
components f^{p) and f~{p), which represent the positive and negative slopes in /, 
respectively, and for which f{p) = f^ijp) — f^{p)- Now let f^{p) and f^{p) be the 
respective fiuence profiles that result when there is gantry motion during delivery. In 
the case of prograde motion, where the effective leaf speed is increased by a positive 
value V, it is shown in Appendix A.l| that for any point p = po 

f^{po-Ap + vf^{po))=f^{po), (7) 

where Ap is an arbitrary constant shift, which could be zero. This means that the 
profile has the same values as the profile, but it is stretched in the direction of p. For 
retrograde motion, where v is negative, the corresponding equation (A. 21) is derived in 



Appendix A. 2 and looks a little more complicated because leaves can move backwards 



in this case. However, because the stretching or compression is generally different for 
and /~, the impact on / = /"*" — /~ can be more substantial. 
Let us now illustrate the findings above for the example case shown in figure [3j The 
dose distribution Dq(x, y) from the first of the 10 beams used injsj^b) (the vertical one at 
00 = 0°), is shown in figure |4|^a). It is composed of the following Chebyshev polynomials: 
DQ{x,y) = (0.935To(x) - 0.018T2(a;) - 0.312T4(a;) + 0.076T6(a;) - 0.262T8(a;))/10. Note 
that there are only polynomials of even degree in this beam profile because of the left- 
right symmetry of the example case. The effect of gantry motion is simulated over 
an arc segment of 18°, from —9° to 9°, during the delivery of this beam. The gantry 
rotates counter-clockwise. The leaf motion is from left to right. Thus, in the upper half 
of figure |4](b), the resulting motion of the field edges in the x-direction is reduced by 
the gantry motion (retrograde motion), and the effect is modeled with equation ( |A.21 ) 
with negative v. In the lower half, which exhibits prograde motion, equation ([T]) with 
positive V is applied. 

The overall result is that, with motion, the characteristics of the fiuence profile are 
preserved, and the effect of motion appears to be similar to using a divergent beam. 
In Single-Arc, the effect is reduced because the MLC moves back and forth during a 
rotation of the gantry, i.e., there is an alternation between retrograde and prograde 
motion. The total dose distribution, which results from combining the 10 arc segments 
according to the 10 beams used in figure [sj^b), is shown in figure [5| Comparing the two 
figures and|5|, one can see that the effects on the dose distribution within the 

unit circle encompassing the target and the critical structure are small. A numerical 
analysis reveals that the average (RMS) difference is 0.02, i.e., 2% of the prescribed 



I This effect is quite different in the case of beam modulator ('compensator') based IMRT and 
tomotherapy, in which gantry motion leads to a 'smearing' of the fiuence distribution. 
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Leaf sweep delivery, gantry arc 0.0 degrees 




X 



(a) 



Leaf sweep delivery, gantry arc 18.0 degrees 




X 



Figure 4: The effect of gantry motion during leaf sweep IMRT delivery is illustrated 
here for one of the 10 beams that produce the dose distribution shown in figure [sj^b). 
Specifically, the vertical beam (0o = 0°) from that example is considered. The dose from 
that beam without motion during delivery is shown in (a) above. Here the beam comes 
straight from the top. In (b) gantry motion is simulated in counter-clockwise direction 
from (j) = —9° to 9° during leaf sweep IMRT delivery of that beam. The leaves move 
from the left to the right. In the upper half of part (b) (for positive y), the resulting 
motion of the field edges in the x-direction is reduced by the gantry motion (retrograde 
motion), whereas in the lower half (negative y) it is increased (prograde motion). At 
y = there is no effect of the gantry motion. 
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dose. Outside of the unit circle, some noticeable smoothing occurs in Single-Arc, and 
the differences are bigger. One can also see that the Single-Arc solution is slightly 
asymmetric. Incidentally, if an odd number of beams had been used in this example, 
the solution would have been perfectly left-right symmetric. 



Single- Arc, 180° 




motion in the Single- Arc mode. The first of the 10 beams is delivered with MLC leaf 
sweep while the gantry rotates from —9° to 9°, the second during the gantry arc segment 
from 9° to 27°, and so forth. The total Single- Arc is being delivered over 180°, from 
-9° to 171°. 



7. Discussion 

In this paper an attempt was made to approach the beam angle problem in IMRT from 
a fundamental theoretical direction. The objective was to derive analytical expressions 
that would allow one to draw some generally applicable, rather than case- specific, 
conclusions. Indeed this objective has been met. However, it could only be accomplished 
at the price of a strongly simplified model. The question arises: what is the validity of 
the findings in the real world? I claim that the results remain valid if one uses more 
reahstic dose models, and provide the following arguments in support of this statement. 
The physical effects that were neglected in the oversimplified model, i.e., dose build-up, 
(exponential) dose fall-off, and beam divergence, lead to smooth variations from the 
model. The real world can therefore be modeled by adding small-degree polynomials to 
the simple model. As long as the degree of these correction polynomials does not exceed 
the maximum degree polynomial used in the model, the conclusions remain valid. 



Number of beams in IMRT, Single-Arc 



13 



With respect to the determination of the required number of beams, I should 
emphasize that equation (|6| can only serve as a rough estimate. First, in the 
interpretation of the Chebyshev polynomials as a series of peaks, I ignored the variation 
of the peak widths (see figure [T| and just considered the average width. A more careful 
analysis would have to take the variation in width into account. Secondly, strictly 
speaking equation (|6| is only valid if one allows negative beam intensities. I address the 
negativity issue by use of the zero-offset idea mentioned above. However, the total zero- 
offset depends on the number of beams. A more detailed analysis would have to take this 
second-order effect, which was also mentioned in section |4] above, into account. Third, 
the number calculated with equation ^ is not a strict practical limit on the number 
of beams. In fact, smaller numbers of beams suffice to obtain approximate solutions. 
Along similar lines, it is clear that the required number of beams may also depend on 
the specific geometry of the target volume and the surrounding critical structures, i.e, 
simple cases require less beams. This paper does not provide a case-specific analysis of 
the number of beams, because a general answer is sought after. 

As far as Single- Arc is concerned, a potential criticism of this work is that the theory 
is only two-dimensional. This would not be a problem as long as the real 3D solution 
can be composed of 2D slices. However, in a 2D model one cannot model collimator 
rotation. Others have found that in Single- Arc a 45° collimator rotation leads to better 
results (Otto 2008) than if the MLC leaves move within the transversal planes, as in this 
paper. It is not completely clear where the advantage of a 45° collimator comes from. 
It has been hypothesized that this advantage is due to the fact that, with the collimator 
at 45°, in parallel opposed beams the leaves of the MLC move in orthogonal directions. 
In other words, parallel opposed beams are not redundant (unlike in the present model 
with a flat dose profile, where they would indeed be redundant). It should be kept 
in mind, however, that the model here assumes a gantry rotation of only about 180°. 
Hence, the delivery could be faster, maybe twice as fast, as a Single- Arc treatment over 
360° with a 45° collimator rotation. The 45° collimator rotation may be of no advantage 
if the Single-Arc angle is restricted to 180°. 

With respect to future developments, the idea to decompose beam profiles into 
well-behaved smooth functions, such as Chebyshev polynomials, may not only be of 
interest for the analysis of the required number of beams, but it could also be beneficial 
for general IMRT treatment planning and optimization. 

8. Conclusions 

In IMRT the required number of beams depends directly on the complexity of the fiuence 
(intensity) profiles that can be delivered within the physical and technical constraints 
of the treatment machine. One can measure the fiuence complexity by the degree of the 
Chebyshev polynomial that is required to approximate the profile. The required number 
of beams then equals the highest degree polynomial plus 1. Rotating the beams over 
short arc segments between the fixed beam positions merely causes the dose distribution 
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to fan out, but does not normally alter it substantially. This, in combination with the 
first conclusion above, serves as a theoretical justification of single-arc IMRT (VMAT). 



Appendix A. The impact of gantry rotation on leaf sweep IMRT delivery 

Let us first review the mathematical basis of leaf sweep MLC dehvery. Consider a (non- 
negative) fluence profile f{p) in the interval [pa^Pb] such that f{pa) — fipb) — 0. One 
can always decompose / into two components, one of which, /+, contains all positive 
slopes and the other, /~, contains all negative slopes: 

f{p)-rip)-rip) (A.i) 



where 



p 

1 






dp 



dp 



Pa 



and 



rip) = 







/- 


dp 



dp. 



(A.2) 



(A.3) 



Here [gY^ equals g \i g > and otherwise. Similarly, [g]" equals if < 0, and 
otherwise. Note that both f~^{p) and /"(p) are monotonically increasing functions. 

Let us first consider f~^{p). For any argument p — Po one can write f^ipo) in the 
form: 

Pb Pb 
/+(po) = / t{p)5{po-p)dp= I f+{p)6{p-po)dp. (A.4) 



Pa Pa 

Using integration by parts, i.e., J Fgdp — FG — J fGdp, one obtains further: 

Pi, 

]Pb>po f df 



r{po)^[r{p)H{p-po)]ii 



dp 



H(p-po)dp, 



(A.5) 



Pa 



because the integral of a delta function is the Heaviside step function, H. With the 

identity H{p — po) — 1 — H{pq — p) one can then easily show that 

Pb 

dt 



fHPo) 



dp 



H{po -p)dp 



(A.6) 



H{po-p+{f))df. 



(A.7) 



The case for / (po) can be made completely analogously (simply replace with 
/~ and p"*" with p~), yielding: 

/max 

/(po)=/+(po)-r(po)= j H{po-p+{f))-H{po-p-{f))df.{A.8) 
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This is indeed the standard uni-directional leaf sequencing algorithm in which the 
trajectory of the left (A) leaf is given by p^{f) and the trajectory of the right (B) 
leaf by p^{f). Here one should identify / with time, t. 



Appendix A.l. Prograde motion 

For the sake of simplicity I will assume that gantry motion adds constant speed, v, to 
the speed of the leaves. Variable speed can be considered but at the expense of a more 
complex notation. Let us also consider a constant shift, Ap. In the following, the tilde 
stands for gantry motion, so: 

P^{f)=P^{f)-Ap + vf. (A.9) 

Consider first the case where the gantry motion is in the direction of the moving MLC, 



that is, in the +p direction, such that v is positive. From equation (A.7) (tilde version) 
one obtains, using integration by parts: 

/max 

f^{po)= J H{po-p^{f))df (A.IO) 



Pb 

= I H {po - pHp)) ^dp (A.ll) 

Pa 

= [H{po-pHp))r{p)]Z 



Pb 



5(Po-p+(p)) i-l-v^] rip) dp. (A.12) 



dp 

pa 

Here p^{p) = p — Ap + vf^{p) and ^—1 — "^^j = The first term in equation 



(A.12) is zero. Now use the relationship 5{g{p)) = J2i \dgJdp{p-)\ ' where pi are the roots 
of g{p). In this case there is one root of Pq — p'^{p), which is p = {p'^)'^{po)- The 



denominator of the relationship for 6{g{p)) cancels out the term — ^—1 — 



m 



equation (A.12). This leaves us with: 

Pb 

r (Po) = [s{p- {p^r'ipo)) /+(p) dp (A.13) 



Pa 



= t {{p^)-\po)) , (A.14) 
which can be written in the form: 

r iP^iPo)) = r {po - Ap + vripo)) = npo). (a.i5) 

Similarly, 

/- {po -Ap + vf-{po)) = ripo). (A.16) 
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This means that is just with variable shifts in the direction of p. However, because 
those shifts can be different for and /~, the impact on / = — /~ can be more 
substantiaL 

An even simpler way to understand the impact of gantry motion is by looking at 



the gradient of the fluence profile. From equation (A. 9) one can see that 
dp^ dp^ 



df df 

or 



+ v (A.17) 



df^ df^ I ( df^\ ,^ , 

^ ''l + v^] (A.18) 



dp dp / \ dp 

Hence, in areas with large gradients of the gradient of will converge to 1/v. 
Areas with small gradients will not be affected. 

Appendix A. 2. Retrograde motion 

Now let us look at the case where the gantry motion is in the opposite direction of 
the moving MLC, that is, in the —p direction, which means that v is negative. This 
is more complicated because now the MLC leaves can move in both directions, i.e.. 



p^if) is no longer necessarily monotonic. As one can see from equation (A. 9), leaves 
move backwards, i.e., p^{f) has a negative slope (and so does, consequently, p^ij))) if 
df^ /dp > —1/v. Let us first look at f~^{p) again. 

From equation ( |A.12 ) it follows that in this case 



r (Po) = fm..H {po - p^ipb)) + Yl (A-19) 

i 

where the Pi are, again, the roots of g{p) = Pq —p'^{p), and Sj = —1 if the slope of p^{p) 
is negative at p = pi, and Si = +1 otherwise. 

Now the task is to calculate f~^{pi). One can keep this problem tractable by 
subdividing the fluence interval [0, /max] into n contiguous sub-intervals [ft^fi+i] such 
that = and /^^^ = /max- The sub-intervals are chosen such that in the i-th sub- 
interval p'^if), and thus p^{p), is either monotonically increasing or decreasing, in an 
alternating sequence. The corresponding p intervals are [pf ,pt^i],i = 1, . . . ,n, where 
pf = Pa and Pn+i = Pb- It is now clear that there can be at most n roots p,. 

Now define 

fi [p- Ap + vr{p)} = < . (A.20) 

I U otherwise, 

which are the monotonic segments with 'zero-padding'. One can see that {p^{p)) = 
f^ip). Hence, at the roots pi one gets f^{pi) = iP^iPi)) = /i^(Po)- This yields the 
final result 

/+(P0) = fra..H {po - p^ipb)) + ■'^^ftiP0)■ (A.21) 
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One can understand this result intuitively as follows: Assume that for any given 
point po the MLC leaf crosses that point pq twice, once, in segment n — 1, in the —p 
direction, and then, in segment n, in the +p direction. The total fluence is then simply 
fniPo) ~ fn-i{Po), because this is how long point po is irradiated. If the leaf crosses the 
point multiple times, the total fluence is exactly as given by equation (A. 21) above. 

In words, equation (A. 21 ) means that the fluence delivered is the sum of compressed 
segments from the original (no motion) fluence profile, and the negative terms therein 
ensure a contiguous transition between those segments. For f^{p) one gets the result 
simply by replacing all '+' in the exponent with '— '. The total overall fluence is 

/M = /+(Po)-rM. 
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